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We study the nature of the superconducting state, the origin of d-wave pairing, and elementary 
excitations of a resonating valence bond (RVB) superconductor. We show that the phase string 
formulation of the t — J model leads to confinement of bare spinon and holon excitations in the 
superconducting state, though the vacuum is described by the RVB state. Nodal quasiparticles are 
i obtained as composite excitations of spinon and holon excitations. The d-wave pairing symmetry is 

shown to arise from short range antiferromagnetic correlations. 

I. INTRODUCTION 

This paper concerns the nature of the superconducting state and the excitation spectrum obtained from the bosonic 
resonating valence bond (b-RVB) theory of the t — J model. The t — J model, arguably the simplest description of a 
^ , doped Mott insulator is given by 
5h ■ 

H = "* E 4(1 - n M )c, a {\ - nj- a ) + j£ [Si • Sj - . (1) 

\ (ij) 

The b-RVB theory of the Hamiltonian (Q) yields a consistent description of the doped antiferromagnetic Mott insu- 
re I lator M. At half filling, the theory reduces to the Schwinger boson representation of the Heisenberg model At 
£5 ■ small doping concentrations, the theory describes a cluster spin glass phase, which eventually makes way to a super- 
conducting phase M. The theory also exhibits a spontaneous vortex phase where the holons are (Bose) condensed 
i^i - and the spinons are perceived as half vortices by the holon field H. Recently, we considered a Ginzburg-Landau 
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formulation of this theory and showed how g£ quantization arises in the superconducting state. A spinon, perceived 
as a half vortex by the holon field, is trapped within a vortex core, thereby leading to a 2e flux quantum In this 
paper, we present a microscopic description of the superconducting state at zero temperature. In particular, we show 
explicitly that when the holons are condensed, it costs a logarithmically divergent energy to create a bare spinon 
excitation. Therefore, a bare spinon excitation is confined. We show that energy considerations only permit neutral 
5 = 1 and quasiparticle excitations. We also demonstrate that the short range antiferromagnetic (AF) correlations 
lead to d-wave pairing symmetry. The results presented in this paper illustrate how a nodal quasiparticle arises from 
an RVB vacuum, at low energies. This paper complements the macroscopic description of the superconducting state, 
presented in our earlier work || . 

The paper is organized as follows. For completeness, we devote the next section to a brief summary of the b- 
RVB formulation of the t — J model ([!]) . This leads to an effective Hamiltonian for the spin and charge degrees of 
freedom which can be solved within a mean field approximation. The basic results used repeatedly in this paper 
are summarized. In the following section, we discuss the nature of the superconducting state, and show how d - 
wave pairing arises in this theory. The excitation spectrum is studied in Section IV. We show that bare spinon (or 
holon) excitations are confined in the superconducting state, and that only neutral S = 1 and composite excitations 
(that carry both spin and charge) are allowed in the superconducting state. The latter would be the RVB analog 
of quasiparticles in a BCS superconductor. In the final section, we summarize our results and discuss them in the 
context of the complete phase diagram obtained within the b-RVB theory. 



II. BOSONIC FORMULATION OF THE T - J MODEL 

The b-RVB formulation can be regarded as a bosonization scheme for the t — J model (Q). The electron operator, 
Ci a is written as 

c la = h\b l(J e^° (2) 

where the holon h\ a and spinon bi a operators represent bosonic fields. The holon and spinon fields satisfy the 
constraint, h\hi + J2a^l<7^ = 1> as m tne fermion or d-RVB theories ||. However, the presence of the nontrivial 
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phase string operator, 6^, in (g), makes this theory very different from the d- RVB theories. The phase string 
operator leads to nonlocal correlations between the charge and spin degrees of freedom on account of hole motion in 
a background with short range AF correlations. The fermion statistics of ci a is also realized through the phase field 

The Klein factor (<j) Nh guarantees anticommutation relations between electron operators with different spins. Unless 
otherwise stated, this factor will be omitted as it does not play any role in the spin, transport, and single particle 
channels. The phase string operator 0^ rmg , is a nonlocal operator given by 

C rinS ^[^-^] , (4) 

where 

" ( " ) , (5) 



and 

$,* = 5>(0n?. (6) 

Here, n\ a and are spinon and holon number operators, respectively. The quantity 0i(l) is defined as 0i(l) = 
Im In (zi — z{), where zi = x\ + iyi is the complex coordinate of a lattice site /. 

Using the decomposition (||), we find that the hopping term H t in ([!]) can be written as 

H t = -tJ2 (e a «"<) h% (f A n) b] a b ia + h.c. , (7) 

where 



(8) 



and 



^• = ~EM-^(0](E™t) > 

= 5 E^W -**(*)] ■ (9) 



The lattice gauge field Afj describes fictitious fluxoids bound to spinons satisfying 

E^ = ±-EKt-<] . ( 10 ) 



ie<- 



for an arbitrary closed path c (assuming no spinons along c) . The sign ± on the right hand side of the above equation 
refers to the clockwise or counter clockwise choices for the path c. The quantity 0?- describes a uniform flux threading 
through the 2D plane with a strength ir per plaquette, J2a 4>% — ^ z7T - The gauge field is given by holons, 



4 = iE^(0-^(0]«?. (11) 
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and describes fictitious fluxoids bound to holons, 

E4 = ±E<- ( 12 ) 



The superexchange term Hj in (Q) is given by 

J 



(ij)aa' 



E (^) ( eiff ' 4i ) . ( 13 ) 



From equations (f?l) and ( |l3| ) , we see that the b-RVB theory can be thought of as a 7r-flux theory where the n fluxoids 
are bound to the constituent particles, i.e., the spinons and the holons. 
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A. Mean field solution 



The Hamiltonians for the charge and spin degrees of freedom, (Q) and ( |l3| ) respectively, can be solved within a 
generalized mean field approximation. We expect the presence of the nonlocal fields Afj and to cause residual 
effects of the short range AF correlations within such a mean field treatment. We summarize below, a few results 
that will be used in this paper. We refer the reader to an earlier paper by one of us for details 0. 

A mean field decoupling leads to an effective Hamiltonian, iJ c g = Hh + H s , where the effective holon Hamiltonian 

H h = -t h £ (e^*-^) h\h 3 + h.c. , (14) 
with th ~ t. The effective spinon Hamiltonian is given by 



H s = -J, (^«) W-„ + h.c + A (j2 b l b ™ 

(ij)<7 \ i<7 



b i<T -(l-S)N\ , (15) 



where J s and J are related by the generalized RVB order parameter A s as J s = %- J. In the above expression, 



A'^/^^/e-^b^bj-J) . (16) 



The single occupancy constraint is enforced on the average by the Lagrange multiplier A in H s . Note that the 
RVB order parameter A s above, incorporates the link field A\j , and is invariant under the following transformation: 
A^j — * A^j +9i — 6j, and b ia — ► b icr e laei . Equations ( |l4| ) and (|iq ) define the effective Hamiltonian, H c ff . All the results 
presented in this paper are obtained using this Hamiltonian. 

The Hamiltonian, Hj can be diagonalized by the Bogoliubov transformation, 

bia = ^ (urnaiihmo- ~ V mcr (i)7m-<r) • ( 17 ) 
m 

The Bogoliubov factors, u ma (i) and v ma (i) can be reexpressed in terms of a single particle wave function w m(T (i) as 
u{v) mcr (i) = u(v) m w mcr (i), and w mrT is determined by the eigen equation 

U™ m „(i) = ~ ^ iaA ^w ma (j) . (18) 

j=nn(i) 

The Hamiltonian, H s is then diagonal in the 7 m representation, H s = Y^ ma - ^mJmalmai where (E m ) m i n ~ SJ, as 
shown in M . The RVB order parameter is determined self consistently as 



A h=J2 e itrA * jW ma(i)w* ma {j)(-U m V m ) 



1 + J2(iL a i^ 



(19) 



In the following, we will always consider a real order parameter Afj . In this case, it can be verified that w ma — w* m _ a 
and that (Ay) = Af^ . Finally, it is useful to note that 

(S i .S j )=-i|A s | 2 <0 (20) 

for nearest neighbors, ij. This is consistent with the notion that the b-RVB order parameter essentially characterizes 
the short range AF correlations. 



III. GROUND STATE: SUPERCONDUCTING STATE WITH D-WAVE SYMMETRY 



The b-RVB ground state is underpinned by the RVB order parameter, A^-. At low temperatures, the holons are 
condensed. At T = 0, the ground state is superconducting. In this section, we will show that the superconducting 
order parameter has d-wave symmetry. This is due to the short range AF correlations. When a hole moves from a 
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lattice site i + x to i + y, the string operator (||) picks up a minus sign on an average due to the short range AF 
correlations. We will show below that this translates into d-wave pairing of electrons. 

The pair operator between nearest neighbor sites ij, A|p = J2a vciaCj-o, will now be reexpressed using the boson 
representation (0). For the purposes of this section, we find it convenient to define the quantities, 



$J = 



i(!)Van 6 



la > 



and 



Then, it is easy to verify that 



where 



$9 



£<mo 



A - 



(_l)J(_l)^ e -i*3-< 



(21) 



(22) 



At low temperature, holons are Bose condensed, viz.. 



hge , where $ is a fixed phase. Then, the 



superconducting order parameter can be written as a product of an amplitude and a phase as 

(M°> = <A2,><e*i(»J+*S)) . 



(23) 



In the Appendix, we prove that the quantity (Ay) is a constant, independent of ij. The important physical quantity 

in this discussion is the phase factor, e' 5 (*<+*j) in (pi]). 

First we note that a real superconducting order parameter will be established only when phase coherence is realized; 
i.e., when spinons are paired up in the ground state such that ^e l5 (* i+ *J')\ ^ 0. Therefore, it should be clear that 
the superconducting transition temperature T c is determined by phase coherence B . 

Besides the phase (de)coherence, the topological phase factor e % 5 ( $i + *j ) in ( pi] ) also determines the d-wave symmetry 
of the order parameter. To see this, let us compare e % ^ ) and e 1 ^ First, we rewrite 

which holds true when two spinons at i and i + rj are annihilated by A° + ~ in ( pl[ ) . We do so because e tAii ^ is well 
defined with (Aj i4 J\ — at T = 0. Therefore, (V^+A » e 2 { (^ii+fi) ) > q in the ground state, and is independent 



of the nn link ii + fj. On the other hand, the difference between e 1 *^ and e l$i +s leads to a phase change. Consider 



_ i[e <+S (i+$)E<r<"'J+s<r- e «+s(^)E -^+eJ+W* 



In the presence of short range AF correlations, one has o^H+ya ~ So- anb i+xa sucn that 

= -1 

since 0i + x{i + 2/ ) — i+ y{i + x) = ±7r. We emphasize that such a such a sign change is generally true on an average in 
an RVB background, and does not require long range AF order. 

Thus, the presence of short range AF correlations among spins causes a d -wave like sign change between ( A^~ 



and / Af^V If one only retains the nn pairing amplitude (Af^J = A°f iii+ fj, with fia+s/fu+y 
momentum space, one has the familiar d-wave form factor 



-1 . then in 



A| C = 2 ( A|C \ (cos k x a - cos k y a) = A° 



(24) 
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IV. ELEMENTARY EXCITATIONS 



In this Section, we consider the confinement of bare spinon and holon excitations. Though the superconducting 
state at T = 0, \^g), is characterized by spinon (RVB) pairing and a holon condensate, it does not necessarily imply 
the presence of free spinon and holon excitations. In fact, we will show that creating a bare spinon or holon excitation 
out of the RVB vacuum costs infinite energy. On the other hand, we will show that neutral S = 1 excitations as well 
as quasiparticle excitations can be created from the ground state with finite energies. 



A. Spin excitations 

A single spinon excitation can be constructed as 7^ l(7 |\f r G). As mentioned earlier, a mean field treatment of 
yields the result, E m ~ 5 J. Therefore, we may be tempted to conclude that the cost of creating a spinon excitation 
is finite. However, the presence of the long range topological gauge field Afj, in the effective Hamiltonian for holons 
( |l5| ) suggests that the creation of a single spinon will be perceived as a half vortex, and consequently would be a 
high energy excitation. To see this, consider the energy cost of such an excitation, (^> G \j ma Hh r yi n(T \'^ G ) — Eq , where 
E' l = (V G \H h \* a ). Now, consider (^ G \j ma e iA ^ ' h\h ]lma ^ G ) . From @, we see that 



I 

Using this, we find 

I 



Noting that — 6j{l) is independent of I when summed over all links < ij > in Hh, and on using the normalization 
condition ^ \ w m<7 {l)\ 2 = 1, we get 



ZJV 

<lj> 

For a nearest neighbor link < ij > far away from any lattice site I, one has 

alsina 



± ]T cos ^(0-0,(0]-! 



i(0-*i(0l 



1 1 J 



ru 



where is the azimuthal angle of the link < ij > and ru ~ is the distance between I and < ij >, and a is the 
lattice constant. Then, a summation over all nearest neighbor links < ij > leads to 

-L £ cos^CO- ^(/)]-l~ln^, 

<ij> 

L being the size of the sample. We thus see that it costs a logarithmically divergent energy to create a single spinon 
excitation. Consequently, single spinon excitations are not allowed in the superconducting ground state at T = 0. We 
reemphasize the crucial role played by the topological gauge field AJ,-, leading to this result. 

The above analysis suggests that a logarithmic divergence in the excitation energy can be avoided by pairing of 
spinons (binding of a vortex-antivortex pair). Such an excitation is nothing but a neutral S — 1 excitation. For 
example, consider a spin excitation created by Sf, Sf\^ G ). At T = 0, the relevant term in Sfl^c) is 7^0-7^/ J^g)- 
The excitation energy can be readily determined from the mean field theory of H s , as 

{^Gh m '. a lm,H sl l al l,_J^ G ) - E° =E m + E m , , 
which is finite. Now let us check the vortex energy in the charge sector, Hh- A straightforward calculation yields 

(*Gf| 7m '_ ff 7maifh7L7^_ ff |*G) - 
<ij> hh 
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In obtaining the above, we have made repeated use of 



hie 2 ^ 



etc. Note that the factor 

~[0<(li)-0i(M]-~[W-0i(la)] 

vanishes when l\ = 1%. It has a dipole configuration when l\ ^ I2, h and I2 being the cores of the spinon 
vortex-antivortex pair. The corresponding energy is no longer logarithmically divergent, since the factor \w ma {li)\ 2 
\w rn ' -<j{h)\ 2 restricts l\ and 1% within a characteristic length scale of the spinon wave packet, l c . Hence, we get (in 
the continuum limit) 



I* 



G) 



^£|«W(MI> m '-«r(fe)| 2 lll 



l r li - r l 2 



hl- 2 



Thus, an S — 1 spin excitation comprising a pair of spinons costs a finite energy, and is an elementary excitation. It 
does not decay into free spinons, due to a logarithmic confining force between the constituent spinons. The excitation 
energy has a lower bound, E g > 2E S , where E s = (£" m )min ~ 6 J. 

To conclude this subsection, we note that an S = 1 excitation created by = (— l)'6||6i|e J * i should also have 
a finite excitation energy, as required by spin rotational symmetry. This can be verified by an explicit calculation. 
In this case, the relevant excitation at T — is of the form lln^lln' i V^g)- Following the same steps as above, 
we find that the logarithmic divergence is avoided by a cancellation between the 2ir vortex and two 7r vortices 
introduced by l m -]l m '+- in the following subsection, we will discuss an analogous effect which leads to a finite energy 
for the quasiparticle excitation. 



B. Quasiparticle excitations 

A quasiparticle (or a quasihole) excitation is obtained as a linear combination of c; CT |<I> G ) and cJ_ |$ G ). We now 
show that the emergence and stability of such an excitation is due to the confinement of bare spinon and holon 
excitations at T = 0. 

From (H), we see that a quasihole excitation has the form 



c iCT |* G )oc^ CT e*[*'- CT *']|tf G ) , 



(25) 



viz., it comprises the holon, spinon and the vortex fields. Let us examine these three components individually. 

First, we have already seen that bi a \^G) would cost an infinite energy in H^. Similarly, a single holon added to the 
ground state will also cost a logarithmically diverging energy in H 8 . It is easily seen that 



(tt G |/nff a fe?|tt G ) 

(*G\hih)\at a ) 



- E s = E° x 



^ £ cos \m)- %(/)]-! 



(26) 



which diverges logarithmically with the lattice size. The third component in (pq) is the vortex field. Consider 



(* G |e+4"*?ff h e-*"*« h |* G ) - E% = E& x 
which is also logarithmically divergent. Similarly, 



i cosi[^(/)-^(0]-l 



G 



(tf G |e-**ifr,e**f|* G ) 



E s n 



2^ E cos |[«*(0- ^(01-1 



which diverges too. Therefore, in the ground state of the b-RVB, the quasiparticle cannot decay into independent 
holon, spinon and vortex excitations as it costs infinite energy. On the other hand, a quasiparticle excitation as a 
composite excitation has a finite energy. In this case, a cancellation of the vortex effects ensues on binding the holon 
and spinon excitations to the vortex fields. 

Consider H s . Clearly, bi a \Q G ) costs a finite energy, and we saw earlier that both /i[|$ G ) and e'*' |4> G ) separately 
cost infinite energies in H s . In contrast, consider the combination hjei^' 1 in (^5|). We will show below that the vortex 
effects cancel. Now, 



H s h}e^'\* G ) = h] 



where the phase twist Xij =1/2 [0i(l) — if (ij)j^ I ancl Xij = ii i = I ot j = I, is responsible for the logarithmic 

energy in (^6|). This is just the statement that the holon is perceived as a half vortex through the term A^. Such a 
vortex effect will be canceled by e5*i , since 



-j a E 1 ' 



bib}.. 



h.c 



Equivalently, 



H„ h}ei*i 



0, such that 



{^ G \h ie -^H s h\ei^\^ G ) = ($ G \hih}H s \<Z> G ) 
is finite. Consequently, Ci a \^ G ) defined in (p5|) must have a finite matrix clement in H s 



,t; 



A similar analysis can be carried out for the matrix elements of Hh- In this case, hl\^ G ) and 1 \*S? G ) cost 
finite energies in H^, but bid^a) and e"^ 0- *' \Jr G ) yield divergent energies. The vortex effects can be cancelled by 
constructing the combination &i cr e - 5 <T *i h ; v { z ^ ff hj b lcr e~%' 7 ' s?l — , which shows that the vortex effect induced by 

bi a in Hh is cancelled by the effect of e~^ a ^ . 

Note that only the vortex phase factor e 5 [* ! ~ CT * i J is single valued. Neither e5*< nor e~5 <T *i are single valued 
by themselves ]l0|| . This means that h[ei^ 1 and bi <J e^^' 7 ^ 1 cannot be independent physical objects. Therefore, the 

holon and spinon fields in ( p5| ) have to be bound to the vortex field e 3 [ $i ] as a whole, such that the quasiparticle 
excitation is well defined and stable in the superconducting state. 

Now that we have demonstrated the stability of the quasiparticle excitation, we determine its energy. Consider 

[if a , C la ] 



-•>> E 

m—nn(i) 



c ia j s e (E & ^-« e ^ L +^ c -) 

m— nn(i) \ ot / 

Js E ^(e^U^^^ 

m—nn(i) \ « 



(m) 



/i.e. 6i, 



+ J S E ^^[E^La 



Each term in the above expression involves three spinon operators and can be linearized using the mean field approx- 
imation ©. Then, 
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[H s , c ia \ ~ 10J s A s c iff +J S ^ /ije^&^e^ (fc^ftt, 

(noting that terms involving e l t e *( m ) vanish after averaging over all possible 0j(m)). The second term on the rhs of 
the above equation can be rewritten, using (pj) and (pl|). We get 

. / aSC\ 



m—nn(i) 



In obtaining the above expression, we have approximated ft] by /i m , since the holons are condensed with an expectation 
value, < h) >= ho (its phase 4> h can be always set as zero due to the gauge choice discussed in Appendix). The 
factor of a in the second term of the rhs comes from commuting the holon field h m with e~ lQm ^" . We have used the 
definition of the superconducting order parameter, 

Af° = ^(/ l IftK i ^^^(-^W m (-l)^eK*H<)) 

CJ 

~ h 2 A s ((- a y(a) m (-l) N »e^ + 0) . 

Then in momentum space, 

[H s , c ka ] « 10A S J s c ka - A^<TC t _ k _ (T , (27) 

where A£ = (^|) . Here, we have used (b l<y b\^ = 1 + (b\ a b la ^ = (3 - <J)/2 and J s = 

On the other hand, the hopping term will generally give rise to a kinetic energy 

[Hh, c kc r] = -(tk - /-i)cka (28) 
with ek ~ th(cosk x + cosfcy). Clearly, a quasiparticle operator can be constructed in the usual manner as 



where u 2 . + v£ = 1. On imposing the condition [iJ s , a k(T ] = E^a^ and using equations (27) and fl28|) , it is 

/ ~ / _ \l/2 _ 

straightforward to get the standard expressions, E 1 ^ = y (et — ^) 2 + (A k ) 2 , Uk = ^7= ( 1 + J , and ^k = 

71 ( 



1/2 

1— %w£ ) sgn(A k ). The first term on the rhs of (27) can be absorbed in the definition of the chemical 



potential, /i. The latter can be determined within mean field theory in the usual manner, by fixing the particle 
number. 

The d-wave gap A k vanishes along the two nodal lines, k x = ±k y . Thus, the quasiparticle spectrum has four nodes 
at ek f — A* = 0. 

It is important to note that even though the confining force between the holon, spinon, and vortex constituents 
within a quasiparticle is weak (logarithmic), the quasiparticle should be a well defined, long lived excitation. Even 
if we neglect the logarithmic energy cost that keeps them bound (which becomes important when the components 
are separated far away), the mean field energies for holon and spinon constituents are zero and E mi respectively. 
Thus, for E^ J < (E m ) . = E s , the quasiparticle energy is always lower than its constituents. This implies that the 
nodal quasiparticle excitation will be observed as a sharp spectral feature. Note that the vortex field in (^5|) also 
provides a confining potential that glues the spin and charge degrees of freedom. This lowers the quasiparticle energy, 
in comparison to the mean field energies of the spinon and holon excitations. 

However, at higher energies, the signature of the composite quasiparticle will be seen in the spectral function. For 
energies uj > E s , the spectral function, A(k, oj) = — ^ImG(k, w)sgn(o;), displays a composite structure in the form of 
a broad peak and spectral weight at higher energies [ |l2||Tc| ] . Therefore, even though low energy quasiparticles in the 
b-RVB superconductor seemingly look like their counterparts in a conventional d-wave BCS superconductor, they are 
easily distinguished from the latter at higher energies, by a broad incoherent structure. This reflects the composite 
nature of the RVB vacuum (i.e., spin charge separation). 

It is clear that the width of the quasiparticle would increase if E^ becomes comparable to, or greater than, E s . 
Typically, this would first happen around the corners of the Brillouin zone, (±7r, 0) and (0, ±7r). In these regions, A k 
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and attain their maximal values. Since the mean field excitation energies of the holon and spinon constituents are 
now comparable to, or smaller than, E^ p , the quasiparticle will be seen as a loosely bound spinon- holon composite. 
This is consistent with results from exact diagonalization |l3|Jl4| ]. 

A similar effect also occurs as the hole concentration is reduced. As 6 — > 0, is finite (recall the factor of in 
the denominator of A£), whereas E s ~ J 6 — > 0. Therefore, at small S, the sharp quasiparticle peak evolves into a 
broad feature quickly as one moves away from the nodal regions to the corners of the Brillouin zone. Eventually, for 
low enough doping, the spectral function becomes very broad away from the nodal points (±7r/2, ±7r/2), and is best 
described in terms of spin charge separation rather than a bound excitation. Elsewhere, we discussed this in detail 
for the case of a single hole doped into the antiferromagnetic Mott insulator pi]. 



V. SUMMARY 



In this paper, we discussed the ground state and the nature of excitations in a bosonic resonating valence bond 
(b-RVB) superconductor. We considered an effective Hamiltonian for the spin and charge degrees of freedom. The 
effective Hamiltonian was derived by rewriting the t — J model using the phase string decomposition. This decom- 
position can be thought of as a bosonization scheme for the t — J model. The ground state at T — is a d-wave 
superconducting state which arises from RVB pairing and holon condensation. 

The physical reason behind d-wave pairing has a simple explanation. Suppose one of the holes forming the Cooper 
pair is located at site i. As the other hole belonging to the pair is transported from site i + x to i + y, say, via i + x + y 
(the shortest path), it picks pick up a phase string, <Ji + x+y<Ji +y , due to an irreparable displacement of the Marshall 
sign that occurs due to hole motion ]lOj . Owing to the short range AF correlations, Oi^x+yVi+y < on the average, 
which implies that the Cooper pair at sites i and i + y changes sign relative to the pair at i and i + x. 

The superconducting condensation in the b-RVB and conventional d-RVB theories ^ both involve the RVB pairing 
of spinons and holon condensation. The main distinction between these two classes of theories lies in the origins of 
kinetic energy frustration. In the b-RVB theory, the kinetic energy of holons is strongly frustrated, due to the short 
range AF correlations that are carefully incorporated. This kinetic energy frustration is realized as the phase string 
effect. A partial recovery of the frustrated kinetic energy occurs at the superconducting transition. However, spinon 
and holon excitations out of the b-RVB vacuum still exhibit the signature of strong mutual frustration between them- 
selves: a single spinon or holon excitation always invites a severe response from the other degrees of freedom, causing 
a logarithmically divergent energy. This is the origin of the low energy quasiparticle, by spin charge recombination. 

We showed that in the superconducting state, S = I excitations composed of spinon pairs have finite energies, 
although single spinons would cost a logarithmic energy. Similarly, a composite quasiparticle excitation, carrying 
both spin and charge, is well defined and has a finite energy. Though the quasiparticle operator we constructed has a 
superficial similarity with the nodal quasiparticle in a conventional d -wave BCS superconductor, there are important 
differences. The first is that neutral S = 1 excitations in the b-RVB state, to leading order, are independent of 
quasiparticle excitations. They occur as bound spinon pairs, whereas in BCS theory, the S = 1 excitation is built 
out of quasiparticles. A second difference is that the composite character of the quasiparticle can be seen in the 
spectral function, at high energies. For low doping, even the low energy quasiparticle is unstable at the Brillouin 
zone corners. It should be noted that the quasiparticle in the b-RVB theory becomes incoherent at T c . As discussed 
in Ref . B , the superconducting transition occurs when spinons are deconfined such that free spinon- vortices destroy 
the phase coherence of the order parameter (]23) . The proliferation of spinon vortices leads to the screening of the 
logarithmic confining potential, and the quasiparticle decays into its constituents. These differences in the excitation 
spectra reflect the fundamental difference between a BCS superconductor and an RVB superconductor. The former 
is constructed from a Fermi sea of electrons, whereas the latter is built from an RVB vacuum which describes a doped 
antiferromagnetic Mott insulator. 

Though spin charge separation is the key feature in the description of the doped Mott insulator, an important point 
we discussed in this paper is that bare holons and spinons need not necessarily exist as excitations of the RVB vacuum. 
The fact that a single holon or spinon cannot be created or annihilated in the system at T = is consistent with 
the requirement of the no double occupancy constraint at the global level. In contrast, an 5 = 1 spin excitation or a 
quasiparticle excitation does not violate such a constraint. The confinement of single spinon and holon excitations, 
along with the emergence of S = 1 spin excitations and quasiparticles as the low lying elementary excitations in 
the b-RVB theory has furthered our understanding of the concepts of the RVB state and spin charge separation. 
In fact, deconfinement of spinons and holons plays a key role in determining the superconducting transition and in 
the description of the high temperature phases of the b-RVB theory. The phase string decomposition (|^) provides a 
natural realization of a two dimensional fermion-boson transmutation envisaged by Baskaran, during the early years 
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of the RVB theory . Recently, Sarker has discussed the recombination of spin and charge degrees of freedom in 
the superconducting state of the t — J model, using a Schwinger boson representation |l7j . Spin-charge recombination 
in the superconducting phase was also conjectured previously by Wen, Lee and others Q within the slave-boson 
framework. In contrast, fractionalization of electrons or spin charge separation of the elementary excitations persists 
in the superconducting state obtained within Z2 gauge theory Jl9| . 

Finally, in this paper, we were mainly interested in the stability and finite energy of the nodal quasiparticle. We 
showed that the composite character of this excitation can be observed clearly as a function of doping. As the doping 
is reduced, the quasiparticle peak around the corners of the Brillouin zone would broaden and the spectral function 
is better described in terms of spin charge separation. A quantitative analysis of this phenomenon in terms of the 
life time of the quasiparticle excitation, and a comparison with experimental results available on the evolution of the 
Fermi surface with doping, will be discussed elsewhere. 
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APPENDIX A: APPENDIX: EVALUATION OF ^A?j, 

The amplitude of the superconducting order parameter in ( p3[ ) can be expressed by 

(A°.)= 5 ^A S (Al) 

with 

g .. = (_i)i(_ 1 )W* e -i(*S+ 2 *?)-<[*?i+(*« h -«)] (A2) 
The link quantity g^ defined in (|A^) may be rewritten as: 

9ij = (-iy(-l) N »e~^+ 2 ^e-^ (A3) 

with % = <P% + - 

By using 9i(j) — 6j(i) — ±7r, one can generally show that 

9ij = goe-^c^e-^i 

where go = (— (— l) Nh e ~ 4( -*Jo +2< ^io-' and c represents an arbitrary path composed of nn links on the lattice which 
connects the site j and a reference site jo- Note that, as ^ n 20^- = ^ n 20^- = ±2tt, g^ is independent of the choice 
of the path c. 

Once jo is fixed, g^ on every link will be determined. For example, on the horizontal bonds connected to jo'- {jo, 
jo ± x), (j , jo ± 2x), (j , jo ± 3x), ... 



(A4) 



9j ,j ±x - 9oe 
g j0±%jo± 2Z = 5o e~ j(2 ^o ± -, 3 o ±S + 2 < ±S .o)-^ o ± ^o ±2 



Similarly on the vertical bonds connected to jo: (jo, jo ± y), (jo, jo ± 2y), (jo, jo ± 3y), 



9j ,jo±y = 9oe l2 ^o±i,™ <*i .io±s = g e i<k io,}o±9 

(A5) 



9j ±y,j a ±2y = 3oe - l(2 ^o±^-o±S+ 2 ^o±i,,o)-<±S- J o±26 
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Now let us consider the choice of the gauge in ^ ■ which satisfies 

= ±7T 

□ 

The holon Hamiltonian in the holon condensed phase reduces to 

H h = -t h hlJ2e\- At ^ +h.c. 

(ij) 

Note that in the ground state all spinons are paired up such that Afj is cancelled out and (e 
which is generally independent of {ij). Then 



{H h ) = -t h hl(e iA '*)Y,e- i 

(ij) 

which is optimized at the staggered- 7r-flux gauge choice of 

±(-1)*tt/4 



'*) + h.c. 



J? 

"TO 

Pii+y 



T(-1)V4 



satisfying (A6). In (A9), there are two choices of signs, ±, for the staggered ir flux at a given plaquette. 
On averaging over these two equivalent gauge choices, we get 



9 jo Jo±3 



9j ±x,j ±2x - 9o \e 



9o 



9o 



(A6) 
(A7) 

(A8) 



(A9) 



In general 



9ij = 9o/V2 



(A10) 
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